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Making use of exact results and quantum Monte Carlo data for the entanglement of formation, 
we show that the ground state of anisotropic two-dimensional S—1/2 antiferromagnets in a uniform 
field takes the classical-like form of a product state for a particular value and orientation of the 
field, at which the purely quantum correlations due to entanglement disappear. Analytical expres- 
sions for the energy and the form of such states are given, and a novel type of exactly solvable 
two-dimensional quantum models is therefore singled out. Moreover, we show that the field-induced 
quantum phase transition present in the models is unambiguously characterized by a cusp minimum 
in the pairwise-to-global entanglement ratio R, marking the quantum-critical enhancement of mul- 
tipartite entanglement. A detailed discussion is provided on the universality of the cusp in ii as a 
signature of quantum critical behavior entirely based on entanglement. 

PACS numbers: 03.67.Mn, 75.10.Jm, 73.43.Nq, 05.30.-d 
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The study of entanglement in quantum many-body 
systems is an emerging field of research which promises to 
shed new light on our understanding of complex quantum 
models. Entanglement represents indeed a unique form 
of correlation that quantum states do not share with their 
classical counterparts, and it thus accounts for the purely 
quantum aspects of the many-body behavior. Given that 
collective phenomena show up in a dramatic fashion at 
phase transitions, the study of entanglement in quantum 
critical systems represents an intriguing subject 
The main focus has been so far the study of spin-1/2 
chains, mainly because they provide paradigmatic ex- 
amples of exactly solvable quantum systems showing a 
quantum phase transition 3]. Only few studies Q have 
explored models in two dimensions and higher, and it is 
therefore hard to tell to which extent the behavior of en- 
tanglement observed at one-dimensional (ID) transitions 
reflects universal critical features (5|. 

A similarly intriguing question in this framework is: 
can we learn anything from entanglement that we did not 
know already from conventional equilibrium obscrvables? 
In general entanglement at T=0 gives a unique insight 
on the global properties of the ground-state wavefunc- 
tion, measuring in a sense how far the quantum ground 
state sits from any possible classical counterpart. There- 
fore entanglement is a powerful tool to detect the occur- 
rence, concrete albeit surprising, of classical-like states in 
strongly interacting quantum systems. 

The aim of this letter is to show that a proper anal- 
ysis of entanglement estimators, applied to interacting 
quantum spin systems in arbitrary dimensions, not only 
exhibits universal features at continuous quantum phase 
transitions, but also unveils the occurrence of unexpected 



ground-state features. 

We consider the two-dimensional (2D) antiferromag- 
netic spin-1/2 XYZ model in a uniform magnetic field: 
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where J>0 is the exchange coupling, (ij) runs over the 
pairs of nearest neighbors, and h=g^,^H/J is the re- 
duced magnetic field. For the sake of simplicity we 
will hereafter understand the canonical transformation 
5r'^^(-l)^5'f with / = 1(2) for i belonging to sublat- 
tice 1(2), which transforms the coupling in the xy plane 
from antiferromagnetic to ferromagnetic. Eq.QJ is the 
most general Hamiltonian for an anisotropic spin-1/2 sys- 
tem with exchange spin-spin interactions. However, as 
real compounds usually display axial symmetry, we will 
henceforth consider either Ay=l or Az=l. Moreover, we 
will apply the field along the z-axis, i.e. ft.=(0, 0, h). The 
case Aj^=l corresponds to the XXZ model in a longitu- 
dinal field, where no genuine quantum phase transition 
occurs upon changing the applied field. 

This paper focuses on the less investigated case Az=l, 
defining the XYX model in a field. Due to the non- 
commutativity of the Zeeman and the exchange term, 
for Ay^l this model is expected to show a field-induced 
quantum phase transition on any D-dimensional bipartite 
lattice, with the universality class of the D-dimcnsional 
Ising model in a transverse field |^. The two cases Aj,<l 
and Ay>l correspond to an easy-plane (EP) and easy- 
axis (EA) behavior, respectively. The ordered phase in 
the EP(EA) case arises by spontaneous simmetry break- 
ing along the x{y) direction, which corresponds to a finite 
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value of the order parameter M^{My) below the critical 
field he- At the transition, long-range correlations are 
destroyed, and the system is left in a partially polarized 
state with field-induced magnetization reaching satura- 
tion only as h~foo. This picture has been verified so far 
in D=l only, both analytically and numerically H,I3- 
We investigate the ground state properties of the 
XYX model on a L x L square lattice making use of 
Stochastic Series Expansion quantum Monte Carlo simu- 
lations using a modified version of the directed- loop algo- 
rithm |l2|,ll0| to account for the low symmetry {Z2 only) 
of the Hamiltonian. With this numerical approach we 
can reach sizes as big as L=28 at very low temperatures 
(T/J=1/2L), such that the T=0 limit of the model, and 
in particular its quantum-critical behavior, is accurately 
captured. 

The main focus of our simulations is on entanglement 
properties, analysed through one-tangle and concurrence. 
The one-tangle quantifies the entanglement between a 
spin and the remainder of the system. It is defined 
as ri=4detp^^'', where p'"^' is the one-site reduced den- 
sity matrix [ill H^ . In terms of the magnetic observ- 
ables, Ti takes the simple form ti=1— 4 ^^(M")'^, where 
a=x,y,z, and M°'={S"). The one-tangle represents a 
global estimate of the entanglement in a translationally 
invariant system, since the vanishing of ri is a neces- 
sary and sufficient condition for the ground state to be 
factorized. The concurrence [l^ quantifies the pairwise 
entanglement between two spins at sites i and j. For 
the model of interest, in absence of spontaneous simme- 
try breaking (SSB), the concurrence takes the form 
Cy =2 max{0,cf-\cff }, where 
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with gf^ = {SfS°'). When SSB occurs, Syljuasen Q has 
shown that Eqs. @-® remain unchanged if the condi- 
tion clj^<C^j^ is satisfied, otherwise Eq. ^ provides an 
upper bound for the actual concurrence. 

One-tangle and concurrence are related by the 
Coffman-Kundu-Wootters (CKW) conjecture 
Ti>T2=J2j^i^ij: which expresses the crucial fact that 
pairwise entanglement does not exhaust the global 
entanglement of the system, as entanglement can also 
be stored in 3-spin correlations, 4-spin correlations, and 
so on. The fact that n-spin entanglement and m-spin 
entanglement with m^n are mutually exclusive is a 
unique feature of entanglement as a form of correlation, 
which puts it at odds with classical correlations. In this 
respect, if the CKW conjecture is verified, the entangle- 
ment ratio R=T2/ti<1 quantifies the relative weight 
of pairwise entanglement, and its deviation from unity 
reveals in turn the relevance of n-spin entanglement 
with n>2. Although indirect, the entanglement ratio is 
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FIG. 1: One-tangle ri and sum of squared concurrences T2 
as a function of the applied field for the 2D 5'= 1/2 XYX 
model with Ay=0.25 and Ay=4. The vanishing of n signals 
the occurrence of an exactly factorized state, while its spike 
signals the quantum critical point. 



the only accessible estimate of multi-spin entanglement 
that we can systematically implement at the moment. 

In Fig. n we plot ti and T2 for the 2D XYX model 
for two values of Aj^ in the EP and EA case. The most 
striking feature is the non-monotonic behavior of both 
quantities as a function of the field. Although the field 
is expected to suppress quantum fluctuations and entan- 
glement only in the extreme h^oo limit, we observe that 
there exists an intermediate non-trivial value hf at which 
entanglement disappears completely, since ti is vanish- 
ing. This signals an exactly factorized ground state in 
the 2D model, completely unknown before. Above the 
factorizing field h{ both ri and T2 have a steep recov- 
ery, that will later be associated with the occurrence of 
the quantum phase transition. Finally, we observe that 
the CKW conjecture is verified for any field value within 
error bars. 

The behavior of the concurrence is found to be quali- 
tatively the same as that observed in D=l and we do 
not report it here; however we underline that also in D=2 
its range stays extremely short and that, for all values of 
Aj, studied, C^j^ and cff cross each other (and vanish 
simultaneously) at the factorizing field. In particular, for 

h>h{, cfj^''>c|J\ so that Eqs. Ull-lO are accurate in the 
most interesting part of the ordered phase, namely above 
the factorizing field and all the way to the critical point. 

The existence of a factorized ground state in the ID 
XYX model has been exactly proven in Ref. Iisl and its 
entanglement signature has been studied in Ref. 0. The 
evidence for a factorized ground state in D=2, coming 
from our QMC simulations, lead us to the 2D generaliza- 
tion of the exact proof of factorization: We were indeed 
able to demonstrate that a factorized ground state ex- 
ists for the most general Hamiltonian Eq. ^ on any 2D 
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FIG. 2: Phase diagram of the 2D XYX model in a field in the 
easy-plane case (Aj,<l, left panel) and in the easy-axis case 
(Ay>l, right panel). The squares are the QMC estimates for 
hf and he- 



FIG. 3: Entanglement ratio 7? = T2/t\ around the quantum 
critical point of the 2D XYX model (A,, = 0.25 and Ay = 4). 



bipartite lattice. The proof will be soon reported else- 
where 16], but we here outline the essential findings. 

For any value of the anisotropies Ay and Az, there 
exists an ellypsoid in field space 



(l-fA,)(l-fA,) (1+A,)(A,+A,) 

hi 



(l+A3)(Aj,+A,) 



= 4 
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such that, when h lies on its surface, the ground state of 
the corresponding model is factorizcd, |5')={^^i IV'i)- 
The single-spin states jV'i) are eigenstates of (n/ • S), 
nj being the local spin orientation on sublattice /. We 
will hereafter indicate with h{ {factorizing field) the 
field satisfying Eq. at h=hf, the reduced energy 
per site is found to be e^— (l-|-Aj^-|-A2)/2. In the 
particular case of A^^l and h^{Q,0, h), the factoriz- 
ing field takes the simple expression hf=2y/2{l+Ay). 
As for the structure of the ground state, taking 
n/=(cos 0/ sinf?/, sin (/)/ sinfl/, cos^/), the analytical ex- 
pressions for (j)i and 9i are available via the solution of 
a system of linear equations. The local spin orientation 
turns out to be different in the EP and EA cases, be- 
ing 4>i=0, (?!)2=7r, 6>i=6'2=cos-i ^(l+Ay )/2 for A,^ <1, 
and (/)i=7r/2, (/)2=-7r/2, 6'i=6'2= cos^^ ^2/(l+A;^) for 
Aj^>l. The numerical and analytical findings for the lo- 
cation of the factorizing field in the 2D XYX model as a 
function of the anisotropy Ay are summarized in Fig. |21 
In particular, the QMC data for h{ are obtained from 
the vanishing point of the nearest neighbor concurrence, 
which gives the best resolution of the location of the fac- 
torizcd state. In the same figure, we also report the line 
of quantum-critical fields h^, extracted through the linear 
scaling of the spin-spin correlation length, ^^^(/i=/ic)~L. 



The critical scaling of the structure factor 
Sxxi'l=(>)'^L'^^'^~^ at the critical field is fully con- 
sistent with the best estimates for the critical exponents 
of the 2D Isingmodel in a transverse field, z=l 7=1.237, 
and i/=0.630 jaHl- We observe that at the Heisenberg 
point Aj,=l the critical field and the factorizing field 
converge to the (non-critical) saturation field h=A. 

It is not clear yet, not even in the ID case, if the oc- 
currence of a factorized ground state be related with that 
of the quantum phase transition in these models. In this 
respect we notice that the change in the analytical ex- 
pression of the concurrence from C-J' to C-^ ^ at h=h[ 
can be related to a qualitative change in the wavefunc- 
tion at the level of the phase coherence between any 
two spins. The factorized ground state could hence repre- 
sent a crucial step towards a global rearrangement of the 
ground state, in view of the quantum phase transition. 

Let us now move to the analysis of the critical behavior 
of entanglement from the point of view of the entangle- 
ment ratio R which, under the validity of the CKW con- 
jecture, provides a new insight in the relative weights of 
multipartite vs bipartite entanglement. Fig. El shows R 
as a function of the field for both the EP (Aj,=0.25) and 
the EA (A=4) case. We observe that a pronounced dip, 
in the form of a cusp, is exhibited at the critical point, 
signaling a quantum-critical enhancement of multi-spin 
entanglement involving n spins with n>2 at the expenses 
of pairwise entanglement. 

The critical features of R in D=2 are surprisingly anal- 
ogous to those exhibited in the ID case by the same class 
of models . Such observation indicates a universal sce- 
nario for entanglement at a quantum phase transition. 
We specialize here to the case of second-order quantum 
phase transitions, assuming that the correct estimator of 
concurrence is provided by Eqs. ©-lO, which implies 
that, if not vanishing, Cij=gff — j+\gfj'~gfj\. As ri and 
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T2 are continuous (and non- vanishing) functions around a 
quantum phase transition, R is continuous as well, and its 
critical behavior is best predicted from its derivative with 
respect to the dimensionless control parameter, hereafter 
indicated as A (X~h in the transition under investiga- 
tion). 

At the quantum phase transition one of the magneti- 
zations, e.g. Af^, will vanish, Af^ ~ (Ac— A)'' for A^A^T. 
Therefore the A-derivative of ri for A<Ac contains the 
term -9a(-A^'^)^~(Ac-A)^^~^, which diverges to +00 if 
/3<l/2 (i.e., away from the mean-field limit). In pres- 
ence of an applied magnetic field along, e.g. the z axis, 
the derivative —d\{M'^)'^ can also be divergent, this time 
to —00. This is indeed the case for the ID transverse 

Ising model ^ where -dxiM"")^ In(A-Ac). On the 

other hand, in the mean-field (D=oo) limit no diver- 
gence of the derivative is present, but only a disconti- 
nuity jump Therefore we expect the singularity in 
the derivative to become weaker when the dimensional- 
ity is increased, as our 2D results suggest when compared 
with the ID results of Ref. |^. This means that, within 
the Ising universality class and below the upper critical 
dimension, dxTi is dominated by the power-law diver- 
gence to -f-oo of ~d\{M^)'^ . This is a fortiori true in the 
case of phase transitions with (3<l/2 in which the only 
"singular" magnetization is the order parameter. For 
A^A+, instead, dxTi is either finite or, again in pres- 
ence of a field along the z axis, it might be dominated by 
the divergence to —00 of —d\{A'F'Y . Therefore ri has a 
discontinuous derivative at Ac, and, being a continuous 
function, it shows a cusp maximum at the critical point. 
This feature witnesses the quantum critical enhancement 
of global entanglement, as estimated through ri. 

On the other hand, T2 is essentially a sum of spin-spin 
correlators gf°' with li— jj limited by the very short range 
of the concurrence |2| • It is fair to assume that T2 is 
dominated by the nearest-neighbor (n.n.) concurrence, 
which contains the most relevant features for this analy- 
sis. In particular, the nearest-neighbor correlators build- 
ing up the n.n. concurrence are also the fundamental 
ingredients of the energy, along with the magnetization 
in presence of a field. At a T^O continuous transition 
the first derivative of the energy w.r.t. A is continuous, 
so that either the derivative of the correlators is not- 
singular, or, if a singularity shows up in d\M^ , it has to 
be compensated by an equal and opposite singularity in 
one of the n.n. correlators. This suggests that d\T2 is 
at most as singular as 9aM^, with the same sign of the 
possible singularity P, Ql- 

Finally, we consider dxR—{d\T2) / Ti—T2{d\Ti) / {txY . 
For A^AjT, if /3<l/2 and under the conditions of d\M^ 
having a weaker singularity than dx{M^)'^, d\R is domi- 
nated by —dxTi with a power-law divergence to —00. On 
the other hand, for A^A^f, 9a i? might be non-singular 
or, at most, as singular as d\M^ , with a divergence to 



-foo. We therefore obtain that R has derivatives of op- 
posite sign when approaching the critical point from left 
and right, and, being R a continuous function, A=Ac can 
only be a cusp minimum. 

In conclusion, through a systematic study of the en- 
tanglement of formation, we have shown that anisotropic 
5=1/2 quantum Heisenberg antiferromagnets on the 
square lattice in an arbitrarily oriented uniform magnetic 
field display an exactly factorized state for a given field 
intensity, depending on the orientation of the field and 
on the anisotropics. The existence of many real com- 
pounds whose magnetic behaviour is described by the 
Hamiltonian Eq. ^ suggests the possibility of an exper- 
imental analysis of our findings. Indeed, the existence 
of a classical-like ground state with flat correlators and 
absence of quantum fiuctuations should lead to striking 
signatures both in the static and dynamical observables. 

Finally, we have also shown that the field-induced 
quantum phase transition occurring in these systems 
is characterized by critical enhancement of multipartite 
entanglement, in the form of a cusp minimum in the 
pairwise-to-global entanglement ratio R. Based on gen- 
eral scaling arguments, we conclude that such a cusp in R 
is a universal entanglement signature of quantum critical 
points for continuous quantum phase transitions below 
the upper critical dimension. 
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